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Abstract: We present a systematic discussion of Braaten and Pisarski's hard thermal loop 
(HTL) effective theory within the framework of the real-time (Schwinger-Keldysh) formalism. 
As is well known, the standard imaginary-time HTL amplitudes for hot gauge theory express 
the polarization of a medium made out of nonabelian charged point-particles; we show that 
the complete real-time HTL theory includes, in addition, a second set of amplitudes which 
account for Gaussian fluctuations in the charge distributions, but nothing else. We give a 
concise set of graphical rules which generate both set of functions, and discuss its relation to 
classical plasma physics. 
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1. Introduction 



It is a theoretically interesting problem to compute the higher-order perturbative corrections 
received by physical observables in the quark-gluon plasma, especially in view of the present 
experimental program at RHIC [1] and of the future LHC heavy ion program. In particular, 
one would like to understand theoretically what the regime of validity of perturbation theory 
is, for various observables of interest, at least within the setup of a locally thermalized plasma. 

Theoretical studies of the thermodynamic pressure, which is a quantity naturally accessi- 
ble to Euclidean space techniques, have revealed a poor convergence of the perturbative series 
unless the strong coupling constant assumes unrealistically small values, ag ~ 0.1. In [2] this 
behavior was attributed to self-interactions of soft, gT scale, gauge fields {g = y/Airas)- Lit- 
tle is presently known about the corrections received by dynamical quantities, in particular 
those which are leading-order sensitive to the gT scale, such as photon production rates [3-6], 
jet energy loss [6,7], heavy quark energy loss [8,9] and transport coefficients (such as shear 
viscosity) [10]. These quantities all are leading-order sensitive to the gT scale, in the sense 
that they would be logarithmically infrared divergent were the screening effects which arise 
at this scale not properly resummed, and are "dynamical" in the sense that they describe 
real-time physics, making their extraction difficult from Euclidean space correlators (and thus 
from lattice data). 

The plasma effects which arise at the gT scale are usually resummed, in a gauge invariant 
way, by means of Braaten and Pisarski's hard thermal loop (HTL) effective theory [11]. This 
effective theory incorporates, to leading order in g, the effects from the scale T on the scale 
gT, in terms of (nonlocal) effective propagators and vertices. These ingredients can be used 
to build a loop expansion, which, so long as only the soft scale gT enters a problem, is an 
expansion into powers of g. Thus, since the above-listed dynamical quantities are leading- 
order sensitive to the gT scale, one expects them to receive potentially large 0{g) corrections 
from soft loops in the HTL theory. 

Until recently, there had been no calculation of the 0{g) corrections received by any 
of these quantities. We believe there is a certain technical advantage in performing such 
calculations directly in Minkowski space, using real-time (Schwinger-Keldysh) techniques; 
however to the best of our knowledge a systematic discussion of the HTL theory within 
this formalism is presently lacking from the literature. In this paper we provide just such a 
discussion. We have recently applied the theory we present here to a next-to-leading order 
calculation of nonrelativistic heavy quark diffusion [12]. 

This paper is organized as follows. In section || we briefly recall the rules of the real-time 
formalism, in particular within the so-called Keldysh ("r/a") basis. In section ^ we show, 
by a power-counting argument, that the HTL effective theory takes on an especially simple 
form in this basis: the only HTL amplitudes (with external gauge bosons) carry at most 
two a Keldysh indices, though arbitrarily many r indices. In section ^ we compute these 
amplitudes, which have a simple and physically transparent form, and give their generating 
functional. 
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A convenient set of effective graphical rules which generate the real-time HTL theory is 
given in section ^; these rules are essentially a graphical realization of the nonabelian Vlasov 
equations (including, as well, Gaussian fluctuations in the particle distribution functions.) 
We discuss the relationships between our results and previous work, and with the classical 
plasma physics of a gas of point-like nonabelian charges, in section ^. Since our method of 
analysis appears to shed little light on the structure of real-time perturbation theory when 
soft fermions are involved, we leave the analysis of fermionic HTLs to future work. 

2. The real-time formalism 

The real-time formalism allows the description of the dynamical evolution of expectation 
values within some initial state or density matrix (as opposed to "in-out" transition ampli- 
tudes). The formalism is characterized by a doubling of the degrees of freedom: in addition 
to the usual ^^(pi" fields which implement forward time evolution, one should path- integrate 
over a second set of fields, "(;/>2", which implement time evolution backward in time to some 
initial time. We work in the so-called Keldysh r/a basis, obtained via the change of basis 
(j)r = ^{4>i + (t'2) and (j)a = (pi ~ 4>2, and let the initial time at which the system's density 
matrix is defined go to —00. For a review we refer the reader to [13]^. 

Here we merely recall the rules of perturbation theory in this context. In thermal equi- 
librium the propagator is a 2 x 2 matrix, which takes the form 

(GH-G.)(l±n(p«))G«\ ^^^j 
Ga ) 

where n(p^) = l/{e^°^'^ ^1) denote the standard Bose-Einstein and Fermi-Dirac distributions, 
for bosons and fermions respectively. For a free scalar field, the retarded propagator Gr would 
be given as^ Gr{P) = —i/{P'^ + — iep^) in Fourier space. The general form ( |2.lD holds 
nonperturbatively: the propagator is completely determined by Gr. The chief reason for 
using the r/a basis in this work is that Bose-Einstein distributions, which play a key role 
when treating soft physics, appear in only one matrix element of the propagator, and are 
therefore most conveniently managed. 

To perform perturbative calculations, one must sum over the r/a assignments for all 
internal legs in Feynman diagrams, subject to the restriction that the vertices carry an odd 
number of a indices. The vertices having one a index coincide with the standard zero- 
temperature ones, and those having three a indices are smaller by a factor 1/4. External 
r and a fields in correlation functions carry distinct physical meaning: since the difference 
field (pa is analogous to an interaction term which would be added to the Hamiltonian, general 
correlators of a and r fields may be understood in terms of the (retarded) nonlinear response 

""^Our r/a fields correspond to the 1/2 fields of the "physical representation" used by these authors. 
^Our metric is ( — h++), and following finite temperature conventions, we capitalize four-vectors but write 
their components as lowercase. 
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Figure 1: Example of a Feynman diagram in the r/a formalism, with (a) explicit marking of r/a 
indices and (b) our graphical notation. The propagators which carry arrows are retarded, and the cut 
propagator is an rr propagator. 

induced by the a fields on some correlator of r fields [13]. Correlators in which the a field has 
the largest time argument vanish, ^ as t — > cxd. 

We graphically represent Feynman diagrams in the r/a formalism by drawing incoming 
arrows on r fields which enter interaction vertices, and outgoing arrows on a fields. With 
this notation, retarded and advanced propagators bear a single arrow, which points in the 
direction of the flow of time. An rr propagator carries two outgoing arrows, which we separate 
by a "cut": we think of the cut as a place where the time flow can start. Sometimes we will 
omit to draw the arrows on these propagators, which should create no confusion. The (tree- 
level) interaction vertices all have an odd number of outgoing arrows. Our graphical notation 
is illustrated in Fig. |^. Incidentally, this particular diagram yields zero, because it contains a 
closed loop of retarded propagators. 

3. Power counting 

We consider (amputated) vertex functions in which all of the external legs are soft gauge 
bosons; by soft we mean P ~ gT and by hard we mean P T, for all components of P. 
We recall that vertex functions having only one Keldysh a index, often called "fully retarded 
functions" , correspond to a direct analytic continuation of the Euclidean vertex functions [14] , 
and that these HTL amplitudes are of parametric size g^T'^~^, where n is the number of 
external legs [11,15]. Amputated vertex functions with no external a index vanish, because 
correlation functions involving only a fields vanish. There are no HTL amplitudes involving 
external ghost fields [11], at least within the classes of covariant and Coulomb gauges, so 
we will not consider diagrams with external ghost lines. We now show that loop amplitudes 
with n external legs, ria of which bearing Keldysh a indices, can only compete with the 
above-mentioned HTL amplitudes if they are of parametric size g^-^^j^^-n ^ Then we show 
that no hard loop can be parametrically larger than gT^~^, implying that vertex functions 
with rifl > 3 are not part of the HTL theory. We also show that the soft contribution to 
one-loop amplitudes behaves like g'^-'^aj'^^-n g^j^^j produces subleading effects relative to the 
HTLs, although vertex functions with na > 3 are soft-dominated, not hard-dominated. The 
key ingredients entering our power-counting are summarized in Table [l|. 
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Vertex functions with more external Keldysh a indices, when they appear within Feyn- 
man diagrams, tend to be suppressed due to the (absence of) Bose-Einstein factors on the 
propagators connecting to them. To see this, we note that when a vertex function carries an 
a index, the Keldysh index on the remote side of the propagator connecting to it must neces- 
sarily be an r index, because of the absence of an aa propagator. There thus automatically 
exists a corresponding diagram in which the a index on the vertex function is replaced with an 
r, and the ar propagator replaced with an rr propagator; since for soft external momenta the 
latter propagator is larger by one Bose-Einstein factor T/p° ~ 1/g, we see that the original 
vertex function will only compete with the latter if it is parametrically larger by one factor 
1/g. By induction this proves our first claim: amputated vertex functions with n external 
legs, Ua of which carrying a indices, will only compete with the HTL amplitudes with one 
external a index if they arc of parametric size g-^-'^aj^'i-n^ 



Ingredient 


Parametric strength 


Soft (bosonic) retarded propagator 


l/g'^T^ 


Soft (bosonic) rr propagator 




Hard, near light-like propagator 


l/gT^ 


Hard three-point vertex 




Soft three-point vertex 


g'T 


d^Q for hard, near light-like Q 


gT' 


d^Q for soft g 


{gTY 



Table 1: Ingredients which enter our power-counting. 



We now show that it is impossible for a hard loop with n soft external (gauge boson) 
legs to be parametrically larger than gT'~"'. Indeed, let us consider a bosonic loop diagram 
(the conclusion being unchanged for fermionic loops), having n three-point vertices and no 
four-point vertices (it can be shown that the latter get suppressed in general.) The dominant 
contribution from the region of hard loop momentum Q ~ T arises when Q is within gT 
of the light cone, in which case each propagator contributes a large ~ l/gT"^ factor. This 
is because, for generic soft external momenta, no two hard propagators can simultaneously 
become closer to the light cone than gT, and parametrically nothing special happens when 
only one propagator becomes arbitrarily close the the light-cone (because of the corresponding 
measure suppression) ; this region maximizes the number of simultaneously large propagators. 
The estimate 1/gT^ for hard, near light-like, propagators is independent on the hard prop- 
agator being retarded or cut^. The restriction of the integration measure d'^Q to the hard, 
near light-like region produces a factor of gT"^, and counting the n three-point interaction 

^Enforcing the mass-shell condition on some hard propagators, when they arc Grr propagators, can force 
linear combinations of the external momenta to be spacelike, but this does not represent a parametric sup- 
pression of the phase space of the soft external momenta. What we are saying is that, for Q = R + P with 

= and P soft, « 2R ■ P, it is fair to treat l/(2i? • Piier") and 6{2R- P) as parametrically equivalent, 
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vertices each as gT, we find that a hard loop with n external legs can behave at most like 
gj'^-n parametrically. Actually a cancellation occurs when na = 1 so these functions behave 
like g^T^~"', but we will see in the next section that no such cancellation occurs when na = 2. 
Together with the last paragraph, this shows that hard loops with Ua > 3, when appearing 
inside full diagrams, have sublcading effects relative to those with na = 1,2. 

What about soft loops? Using the rules of the r/a formalism, it can be shown that loop 
diagrams with Ua external a indices can contain up to Ha internal cut (rr) propagators. For 
bosonic loops, the presence of Ua Bose-Einstein functions suggests that diagrams with Ua large 
should be soft-dominated, not hard-dominated. This indeed happens: in the soft region, a 
(bosonic) retarded propagator should be estimated as parametrically a cut propagator 

estimated as 1/g'^T'^, a three-point interaction vertex treated as g^T, and the integration 
measure treated as {gT^. Adding up, we find that the soft loop contribution to a bosonic 
loop diagram is parametrically y4-no2i4-m gj^(,g ggf^ loops should actually be evaluated using 
effective HTL-resummed propagators and vertices, which are quite complicated expressions, 
we do not expect parametric cancellations to occur. 

For Ua = 1,2 the soft contribution is down by exactly one power of g relative to the hard 
contribution: these diagrams are truly hard-dominated and deserve to be called "hard thermal 
loops." Hard loops with arbitrary external r/a indices were considered in [16], in which it 
was found that due to a cancellation hard loops with ria = 3 were of parametric order g'^T^-'^ 
(however these authors did not recognize that hard loops with na > 3 only had subleading 
importance in actual calculations). Thus we see that hard loops with ria > S not only have 
subleading effects, they arc also incorrect: the corresponding diagrams are so/t-dominated, 
not hard-dominated. However this poses no problem to power-counting: when present inside 
diagrams, soft loops having arbitrary numbers of external a indices all contribute at the same 
order as the soft contribution to the Ua = 1 loops, e.g. are down by 0{g) relative to the 
HTL contribution. Essentially what happens is that by modifying the external r/a indices 
on a soft loop, one is merely transferring the Bose-Einstein factors back and forth between 
the propagators outside and inside the loop. 

The HTL amplitudes account for the dominant effects from hard particles on soft gauge 
fields, and together with the tree interactions vertices (which are of the same order), they 
can be used to set up an effective theory which contains only the gT scale. We believe that 
the loop expansion within this theory, in the real-time formalism, proceeds in a way entirely 
similar to that in the imaginary-time formalism [11]: as long as only the scale gT contributes, 
each additional loop is suppressed by one power of (7 ^. In general, in a theory with only one 
physical scale, one would expect the loop expansion to be an expansion into powers of y^; in 
the imaginary-time formalism one gets an expansion into powers of g, because each additional 
soft loop can introduce one (and only one) additional Bose-Einstein factor. Within the real- 

''Of course, in practice such an expansion is not expected to hold up to arbitrarily high order, since other 
scales should eventually enter (associated e.g. with some mean free path, or with nonperturbative, infrared 
phj'sics). However, the appearance of a new physical scale would be signaled by a divergence in the HTL 
effective theory, because it contains only one scale. 
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Figure 2: The most general one-loop diagram with one external a index (the outgoing arrow), but 
many r indices: all such diagrams only have one cut propagator. The arrow flow corresponds to time 
flow; we do not specify the nature of the particle running in the loop. 

time formalism, this well-known statement must be modified to the claim that each additional 
loop introduces one (and only one) additional Bose-Einstein factor or HTL vertex functions 
with two external a indices. The latter vertices may be regarded as secretly containing one 
Bose-Einstein factor, since they are larger by a factor 1/g. We will not embark here into a 
general proof of this claim, which would be easy to give using the effective Feynman rules we 
present in section ^; we simply remark that it would be quite surprising for a systematic loop 
expansion to hold within one formalism, and not within another. 

4. Calculation of the hard thermal loops 
4.1 HTLs with one a index 

The gluonic hard thermal loops having only one external Keldysh a index, often called "fully 
retarded functions", can be obtained via a direct analytic continuation of the well-known 
Euclidean ones [11] [15], and thus do not need to be independently recomputed within the 
real-time formalism. Nevertheless, we think it is instructive to briefly describe the Feynman 
diagrams which contribute to them, and how, following [19, 20] they can be evaluated using 
a simple kinetic theory of point-like particles. 

The most general one-loop Feynman diagram with one external a index, but arbitrarily 
many r indices, is illustrated in Fig. |2[ All of these diagrams contain only one cut (rr) prop- 
agator, which occurs at the smallest time in the diagram. In terms of our graphical notation, 
this structure can be tracked down to the fact that inserting an external r field (incoming 
arrow) onto a diagram (involving only bare interaction vertices) , be it onto an existing vertex 
or onto a propagator, never modifies its arrow flow. Physically, this structure means that 
these functions represent the nonlinear response of the expectation value of the gauge current 
(at the a leg) due to a background field (the r legs): the Feynman diagrams mimic the ob- 
vious quantum-mechanical procedure of starting with an initial one-particle density matrix 
(here, the "cut"), evolving it under a background field, and taking the expectation value of 
the current at late times. 

Li [19] the problem of calculating the induced current in a soft background gauge field 
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(mean field) was considered, and reduced to kinetic theory: 

JZiX) = Y.9 [^v^n^iX,p), (4.1) 
DOF ^ ^ 

V ■Dn''{X,p) = gTi'r [t^tj] dv'E'\X)^^^, (4.2) 

wliere = (l,p/p) represents tlie four-velocity of a hard particle and n denote the standard 
Bose-Einstein or Fermi-Dirac distributions. The second equation, to be solved with retarded 
boundary conditions, determines the color-adjoint disturbance n". The concept of point- 
like particles originates from the separation of scales gT <^ T, between the momenta of the 
external gauge field and that of the typical particles which contribute to the induced current: 
the hard particles feel the external field as if they were point particles. The degree of freedom 
count in ( [d.!] ) is: two bosonic degrees of freedom for the gauge/ghost system, four fermionic 
degrees of freedom for Dirac fermions, and two bosonic degrees of freedom for complex scalar 
fields; the form of the resulting equations is the same for all of these particles, and is gauge- 
fixing independent [11]. A collision term is not included in ( |4.2D because collisions are only 
relevant over 1/g'^T time scales; for a more ample discussion we refer to the review [21]. 

Solving for the induced current ([4.1|)-( |4l^ ) with retarded boundary conditions yields the 
term in the generating functional (effective action) of real-time amplitudes which is linear in 
the Keldysh Aa field: 

'dJ^J d*X r dTV^Ar{X)U-"'(X,X -VT)lAAv'E"'lAr]{X -vt), (4.3) 



where U stands for an adjoint Wilson line along the hard particle trajectories. Here we 



have explicitly performed the radial integration in (4.1), leaving only the integration over the 
angle v of the hard particles; in a generic Yang-Mill theory with Np Dirac fermions and Ns 
complex scalars, the degree of freedom counts add up to: 

ml = ^ [Ca + NfTf + NsTs] , (4.4) 

where Ca = -^c and Tp = T5 = ^ in SU(A'^c) gauge theory with matter in the fundamental 
representation. 



We find (4.3) rather physically transparent compared to its Euclidean counterpart [17] 
[18], however compact the latter might be. The equivalence between the analytic continuation 
of the vertex functions derived from ( |4.3| ) and the standard Euclidean ones can be verified 
from the explicit expressions for the induced current (^rEuciidean/'^^ given in [18]. It is a rather 
nontrivial, although necessary, property of these functions that they become symmetrical in all 
of their arguments (including the a leg) when the boundary conditions are made symmetrical, 
by analytically continuing the momenta to imaginary frequencies. However, it is this very 
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asymmetry between the a leg ("induced current") and the r legs ("external fields") of the 
fully retarded vertex functions, at physical values of the momenta, which makes the generating 



functional (O) so simple. We describe the Fourier space amplitudes derived from ( [4.3| ) in 



more detail, in section |^. 

4.2 HTLs with two a indices 

We now compute the hard thermal loops with two external a indices, beginning with the aa 
self-energy. The relevant Feynman diagrams are shown in Fig. |3| (a)-(c); the propagators in 
these diagrams are most conveniently added together under the integration sign: 

Grr{Q)Grr{R) + ^Gr{Q)Gr{R) + ^Ga{Q)Ga{R) ^ \ [G>(Q)G>(/?) + G<{Q)G<{R)\ , 

(4.5) 

the equality holding up to analytic terms which integrate to zero such as G ji{Cli)G a{R) 
("closed loops of retarded propagators"), which we have subtracted in passing to the right- 
hand side. Here the propagators G'*'^ denote the Wightman (unordered) two-point functions, 
and we have employed the identities Gr — Ga = G^ — G^, Grr = hi^^ + G"^)- Equation 



(4.5) states that the aa self-energy is the average of the two Wightman self-energies, a cutting 
pattern which we graphically represent with two parallel lines, as in Fig. |3| (d), in analogy 
with our notation for the Grr propagator. Evaluating the loop for a complex scalar field, for 
definiteness, yields: 

' °^ ' G>{Q)G>{R) + G<{Q)G<{R) 

^ -g" J2 '^r [trtr] [ ^^^^^i^ ' ^) ^M^- ^B(g)(l + ^^(g)) , (4.6) 
DOF ^ 

where on the second line we have used P <^ Q since the integral is saturated for Q ^ T. We 
see that the calculation of this self-energy is relatively simple, compared to the corresponding 
retarded HTL self-energy: since no 0{g) cancellation occurs, all P <C Q approximations 
can be applied directly. Particles of different spins yield similar contributions (up to 0{g) 
corrections) , and one gets the same degree of freedom count as in (|4.1| ) . The self-energy ( [4.6|) 
could also be obtained by means of the KMS (fluctuation-dissipation) relation, which relates 
the aa HTL self-energy to T /p^ times the discontinuity of the retarded HTL self-energy. 

Equation (|4.6|) has a rather straightforward physical interpretation: the double cut in 
diagram (d) describes thermodynamical equal-time fluctuations, 

(n'^(x,q)n^(x',q')) =Tr, M - ^){27,f S^c^- c{)nB{q){l + nB{q)) , (4.7) 



affecting the distribution functions which appear in ( [4.11 ). Expression ( |4.6| ) is precisely (i)^ 
times the Fourier transform of the unequal time correlator derived from (14.71) , by means 
of ballistic propagation v ■ 9n"(X, q) = 0. Considerations of gauge invariance immediately 
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(a) (b) (c) (d) 



Figure 3: (a)-(c): Diagrams contributing to the aa gluon self-energy, with momentum assignments 
displayed on the first diagram, (d) The sum of the three preceding diagrams. The double cut represents 
the average of the Wightman cuts in the two direction. 




Figure 4: The most general form of diagram contributing to gluon HTLs with two a indices but many 
r indices. Diagrams with four point vertices are subleading. The double-cut has the same meaning as 
in Fig. I 

suggest that to obtain the fluctuation functions with external r gauge fields, one should merely 
substitute this free ballistic propagation with a gauge-covariant one, v ■ D[Ar]n"' {X , ci) = 0. 

Not surprisingly, this expectation is borne out by explicit calculations. This is especially 
obvious to see if one begins with a reorganization of the r/a structure of the relevant diagrams, 
in a way analogous to ( |4.5| ) above. Indeed, although a direct application of the rules of the 
Schwinger-Keldysh formalism would yield diagrams which correspond to inserting the external 
r fields onto the diagrams (a)-(c) of Fig. ^, it can be proved that their sum is equivalent to what 
is obtained by inserting the r fields directly onto the simpler diagram (d): all of the diagrams 
thus obtained contain two cut (Wightman) propagators, which appear at the smallest time, 
as illustrated in Fig. ^. Making all small P approximations to the propagators and vertices 
in these diagrams amounts to taking the propagation of the hard particles to be eikonal, e.g. 
given by Wilson lines along their classical trajectories, justifying the procedure mentioned in 
the preceding paragraph of replacing v ■ d with v ■ D[Ar]. 

Thus, upon performing the radial integration in ( [4.6| ), we obtain the generating functional 
for vertex functions with two a indices (with mj-) as in ( |4.4| )): 

r(2) = !^ I^Jd^X 1^ dTv,AriX)U''\X,X -vT)[Ar]v,A:\X -vt) (4.8) 

Although it should be rather obvious, from their manifest gauge-covariance, that the 
vertex functions obtained from (|4.3| ) and (|4.8|) obey a large class of Ward identities, we 
remark that the full HTL effective action F^^^ -1- F^^^ is not strictly invariant under the whole 
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set of Schwinger-Keldysh gauge transformations . Such a strict invariance would require the 
inclusion of terms which are 0{A^) and higher order in the Keldysh fields, which, however, 
we prefer not to include since we argued in section ^ that such terms are not part of the HTL 
effective theory. 



5. Feynman rules for kinetic theory 

We now present a concise set of Feynman rules which generates the complete (gluonic) HTL 
effective theory in the real-time formalism. To obtain these rules we first rewrite the induced 
current 6T/5Aa from (|4.3D into the form: 



(5.1) 



where we have decomposed the electric field as = d^v ■ A — v ■ DA'^ and used J ^^v^ = 5q. 
Expressions of the form {1/v ■ D)S should be understood as the solution of v ■ D = S with 
retarded boundary conditions, e.g. the adjoint Wilson line in 



(a) ^'KSm^SlSU'' =imlS^6^ (d) « < ^ < c =-v^f^bc 

(b) /^vilM^^ =iTv'' (e) 



P ^ -I 

~ V P-^ 



(c) ^^IMyA^ =ip^v^' (f) =27t6{vP) 

Figure 5: Effective Feynman rules for the HTL theory in the r/a formalism. The arrows follow the 
graphical notation for r/a diagrams introduced in section |2[ All two-point functions are proportional 
to the identity in color space, 6°"^, not explicitly shown. A factor {m?jj/T) J must be given to every 
disjoint double line appearing in a diagram. 

The graphical rules given in Fig. ^ reproduce this induced current, which gives the am- 
plitudes with only one external a index (outgoing arrow). Specifically, the A^6q term in (5.1), 



which gives rise to a "contact term" in the retarded self-energy, is mapped to the component 
(a) of Fig. ^. The second, non-local term in ( ^.1] ) is mapped to a class of diagrams, in which 
an incoming gluon first generates a disturbance in the distribution function associated with 
the four- velocity via vertex (c), which is then evolved using the eikonalized retarded prop- 
agator (e) and interaction vertex (d). The current associated with this disturbance sources 
a gauge field via vertex (b). Finally, one must perform the integration over the disturbed 
particle's momentum by adjoining a factor m^j J to all double lines in a diagram. The 
cut eikonal propagator (f), which plays no role in the calculation of the HTL amplitudes with 



^More precisely it is invariant under those infinitesimal gauge transformations whose parameter is a Keldysh 
r field, but not under those for which it is an a field. 
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only one external a index, enters the calculation of amplitudes with two external a indices 
(in which it appears exactly once). By (v ■ P~) we mean 



1 



1 



(5.2) 



vP- 



V ■ P — ie 



The effective propagators and interaction vertices of Fig. ^ 



are to be used in building 



Feynman diagrams according to the standard rules of the r/a formalism; the objects (a)-(c) 
correspond to ar interaction vertices, the propagator (e) is a retarded (ra) propagator and 
(f) represents an rr propagator. The interaction vertex (d) carries arr indices and is the only 
such interaction vertex in this theory; there is no three-point vertex involving an a gluon. 
In practice the self-energies on soft gluon lines must be resummed; when this resummation 
is performed the insertion (a) should be ignored, as well as all diagrams in which only two 
gluons connect to a double line. We have not explicitly shown the tree interaction vertices 
between soft gluons, although these are of the same order as the HTL ones and must be 
included in the effective theory. 

The double lines in our graphical rules may be thought of as two-particle states (alter- 
natively, one-particle density matrices): in general, by "opening up" these double lines the 
one-loop diagrams considered in section Q are recovered, or more precisely, specific sums of 
these diagrams. As we comment on below, these graphical rules represent purely classical 
plasma physics. 

5.1 An example 



Figure 6: EfTective Feynman diagrams contributing to HTLs with three external legs, up to permu- 
tations, (a) HTLs with one a index. HTL. (b) HTLs with two a indices. 

As an example of the application of these rules, we evaluate the three-point HTLs. The 
relevant diagrams are shown in Fig. ^, and upon adding their permutations, the application 
of the rules produces: 




Q, cr, c 



(a) 



(b) 



r, 




(5.3) 



(5.4) 



Our notation for F, which corresponds to (— z) times the amputated Feynman diagrams 
themselves, is that the momenta written before the semicolon are outgoing and the others are 
incoming. In particular, momentum conservation implies P = Q + R and Pi + P2 = Q (see 



Fig. g). Our result for the retarded three-point function agrees with the standard one; we 



have also verified the four-point function. Our result (5.4) for the vertex function with two 
a indices agrees with that of [16], Eq. (Ill.lSe) (one must take into account that their F/j/j^ 
function is 2i times our Taar, that their momenta are incoming, and their metric is different 
from ours.) 

6. Discussion 

In this paper we have given the generating functional for all hard thermal loops in the real- 
time formalism, ( [4.3| ) and ( [4.8] ). The form of these amplitudes is deceptively simple: a classical 
plasma physicist, instructed of the fact that nonabelian charges tend to precess in the presence 
of a gauge field, would have known enough to simply write them down decades ago. Indeed, 
we find that there are only two kinds of HTL amplitudes: amplitudes with only one external 
Keldysh a index, which describe the (nonlinear) polarizability of a medium of nonabelian 
point-charges, and amplitudes with two external Keldysh a indices, which represent current- 
current correlations in such a medium. What we think is most interesting about our findings, 
is the fact that these functions form the complete real-time HTL theory. The essential reason 
for this, discussed in section ^, is that the Keldysh a ("difference") fields, compared to the 
Keldysh r ("average") fields, are unable to take advantage of the large occupation numbers of 
the soft gauge fields, hence diagrams containing more Keldysh a indices naturally tend to be 
subleading. Incidentally, loop amplitudes with more than two external a indices are actually 
soft-dominated, not hard-dominated. 

We have given simple effective Feynman rules, in section ^, which generate the hard 
thermal loop amplitudes; physically these rules are nothing but a graphical representation of 
the nonabelian Vlasov equations (see e.g. [22] or the review [21]). Diagrams not involving 
the cut propagator (f) of Fig. ^, account for the classically induced current ( |5.1| ) due to 
a background mean field, and whenever this (f) propagator appears, its role is to account 



for the Gaussian fiuctuations (4^) of the particle distribution functions (corresponding to 
fluctuations of the ^W" fields in the language of [21]). The importance of such Gaussian 
fluctuations was discussed previously in the context of the hot electroweak theory (see [23]); 
the analysis of the HTL theory given in [24], of which we became aware after this work was 
completed, bears much similarity with ours. 

Our findings extend in a straightforward manner to nonequilibrium setups: one should 



simply replace the distribution functions in (4T) and (4/7) by their time-dependent expres- 



sions. As was argued for in great detail in [25], this procedure will be correct provided the 
naive criterion for it to make sense is satisfied: the distribution functions of the particles 
should be slowly varying on the length and time scale set by the "soft" gauge fields (those 
for which the HTL effects become important). Also, the perturbation theory employed here 
should make sense: this requires a separation of scales between the "soft" scale and the 
momenta of the particles which dominate the integrals (^]^) and (|4.7| ) (in this paper this 
separation is gT <C T). 
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Even though the close relationship between the HTL theory and classical plasma physics 
has been widely recognized for a long time, it was not clear, at least to the author, how 
this understanding could be exploited in the context of the next-to-leading order calculation 
of dynamical quantities like the ones listed in the introduction. The difficulty is that at 
leading order one has to deal not only with soft, classical physics, but also with some truly 
quantum physics. This is the context in which we believe a systematic, a priori fully quantum 
mechanical approach, starting from the Feynman diagrams of the real-time formalism, can 
be most useful. 
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